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We propose a technique to design finite-length irregular low -density parity-check (LDPC) codes over the 
binary-input additive white Gaussian noise (AWGN) channel with good performance in both the waterfall and the 
error floor region. The design process starts from a protograph which embodies a desirable degree distribution. 
I This protograph is then lifted cyclically to a certain block length of interest. The lift is designed carefully to 

. satisfy a certain approximate cycle extrinsic message degree (ACE) spectrum. The target ACE spectrum is one 

with extremal properties, implying a good error floor performance for the designed code. The proposed construction 
results in quasi-cyclic codes which are attractive in practice due to simple encoder and decoder implementation. 
Simulation results are provided to demonstrate the effectiveness of the proposed construction in comparison with 
similar existing constructions. 
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I. INTRODUCTION 

Design of finite-length low-density parity-check (LDPC) codes which perform well in both the wa- 
terfall and the error floor regions is a challenging task. Irregular degree distributions which are op- 
timized to render superior waterfall performance will often result in a high error floor in randomly 
constructed codes. To improve the error floor performance of irregular codes, different approaches have 
been examined. In one direction, iterative decoding algorithms are manipulated to perform better in 
the error floor region, see, e.g., [4]. In another direction, new constructions of LDPC codes are in- 
troduced O, flu, [[I5l, [[61, [ini, [El, [|71, [HI. In this paper, we are interested in the latter approach. 
Our work is closely related to the ideas and constructions introduced in [[T2l . [[T3l . and [[T4l . on one hand, 
and those in |[T5l , |[6l , and [|T||, on the other hand. 

The error floor performance of an LDPC code over the additive while Gaussian noise (AWGN) channel 
is closely tied to graphical objects, referred to as trapping sets IfTOI . A full characterization of dominant 
trapping sets over the AWGN channel, particularly for irregular codes, is not available. It is however 
known that enumerating such sets, in general, is a formidable task [[U, [0. Indirect measures of the 
error floor performance, which are computationally more efficient, have thus been used in the design 
and the analysis of the LDPC codes. It is well-known that cycles in the Tanner graph of the code are 
responsible for the suboptimal performance of iterative decoding algorithms. Different metrics have thus 
been introduced to measure the harmful effect of cycles. The simplest such metric is the length of the 
shortest cycles in the graph, called girth. Using this metric, Hu et al. proposed the progressive edge-growth 
(PEG) algorithm [l5J, which aims at maximizing the girth of the code's Tanner graph. It was however 
observed that not all the cycles of the same length are as harmful in iterative decoding HI 21 . In [1T2| . the 
approximate cycle extrinsic message degree (ACE) of a cycle was introduced as a metric to evaluate the 
harmfulness of a cycle. The larger the ACE, the less harmful a cycle would be among a set of cycles of 
the same length. ACE constrained LDPC codes were also designed in [[T2l which outperformed random 
codes in the error floor region. The ACE metric for a cycle was recently generalized to the ACE spectrum 
of a Tanner graph in [[T4l . where the authors also devised generalized ACE constrained LDPC codes, 
further improving the error floor of the codes. Finally, in [iOl, PEG construction and generalized ACE 
constrained design were combined for even superior error floor performance. 

In 11511 . edge swapping was proposed as a technique to increase the stopping distance of an LDPC 
code, and thus to improve its error floor performance over the binary erasure channel (BEC). Random 
cyclic liftings was also studied in [[TSl and shown to improve the average performance of the ensemble 
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in the error floor region compared to the base code. Ivkovic et al. [[6l applied the same technique of 
edge swapping between two copies of a base LDPC code to eliminate the dominant trapping sets of the 
base code over the binary symmetric channel (BSC). This was then generalized in [1] to cyclic liftings 
of higher degree, where the liftings were designed to eliminate dominant trapping sets of the base code 
by removing the short cycles which are part of the trapping sets. 

In this work, starting from a protograph, we use cyclic liftings to construct LDPC codes. The resulting 
codes are quasi-cyclic and thus attractive from implementation point of view. Cyclic liftings are carefully 
designed to achieve a target ACE spectrum with extremal properties, and thus the name ACE constrained 
cyclic edge swapping for the technique. The extremal ACE spectrum would imply a good error floor 
performance for the code. 

Compared to the constructions of [|T2l . [fT3l . and [fT4l . which do not have any particular structure, the 
proposed construction is quasi-cyclic and thus more implementation friendly. Moreover, the approach to 
achieve a target ACE spectrum is different. In particular, in the proposed construction, we first find a set 
of most vulnerable subgraphs of the protograph whose inverse image in the lifted graph can potentially 
be small cycles with low ACE values. We then carefully assign cyclic permutations to selected edges of 
these subgraphs such that their inverse image satisfies the target ACE spectrum for the lifted graph. 

In comparison with [61 and [HJ, which start from rather large base codes (graphs), we begin with a 
rather small base code. This implies a more compact description for the code which can in turn result 
in simpler encoder and decoder implementation. As a consequence of a smaller base code, to achieve a 
given block length for the lifted code, the lifting degrees in this work are larger compared to those in [|6] 
and [111. Another difference is in the approach to design the cyclic liftings. In [6] and [1], the liftings 
are designed to eliminate dominant trapping sets of the base code, implicitly assuming that such trapping 
sets are known and available. While this assumption may be valid for hard-decision algorithms over the 
BSC, for the soft-decision algorithms over the AWGN channel, the knowledge of dominant trapping sets, 
particularly for irregular codes, is much harder to attain if not infeasible. Moreover, the approach of [|6l 
and [[II, can work only for the cases where the base graph is large enough, and thus sparse enough, to 
allow for the existence of a meaningful set of dominant trapping setsQ This may not be the case for the 
small base graphs which are the subject of this work. 

'To explain this, consider the (155, 64) Tanner code fTTl, which has been used as a base code in constructions of both (5) and (T). This 
code itself is a cyclic lifting of a complete bipartite graph K^,,^ with 5 variable nodes of degree 3 and 3 check nodes of degree 5. While 
one can meaningfully find the dominant trapping sets of the (155, 64) Tanner code for a certain decoding algorithm over a certain channel, 
the same is not possible for the much smaller K^,,^ graph. 
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The remainder of this paper is organized as follows. Section HI] briefly reviews the definitions, notations 
and concepts required in the rest of the paper. The ACE constrained cyclic edge swapping method is 
presented in Section Hill followed by some simulation results in Section |IVl The paper is then concluded 
in Section |Vl 

II. PRELIMINARIES: LDPC CODES, TANNER GRAPHS, CYCLIC LIFTINGS, AND ACE 

SPECTRUM 

A. LDPC Codes and Tanner Graphs 

Consider a binary LDPC code C represented by a Tanner graph G = {VbUVc,E), where H = 
{61, ... , bn} and Vc = {ci, . . . , Cm} are the sets of variable nodes and check nodes, respectively, and 
E is the set of edges. Corresponding to G, we have an m x n parity-check matrix H = [hij] of C, where 
hij = 1 if and only if (iff) the node Cj G Vc is connected to the node bj E Vb in G; or equivalently, iff 
{bj, Ci} G E. If all the nodes in the set Vb have the same degree d,^, and all the nodes in the set Vc have the 
same degree dc, the corresponding LDPC code is called {d^, dc)-regular. Otherwise, it is called irregular. 
For an irregular LDPC code, the degree distributions of variable nodes and check nodes are often described 
by the two polynomials, A(x) = '^fj^2 ^i^^^^ ^iid = J2f=2 Pi^^^^ ^ respectively, where Dy and Dc are 
the maximum variable node degree and the maximum check node degree, respectively, and Aj and pi are 
the percentage of the edges connected to the variable nodes and check nodes of degree i, respectively. 
Alternatively, vectors A = (A2, . . . , Ad„) and p = (p2, • • • , PdJ can be used for the description of the 
degree distributions. In this case, the code is referred to as a (A, p)-irregular code. 

A (directed) walk of length n in a graph G is a non-empty alternating sequence f oCif 1 . . . Vn-iCnVn of 
nodes and edges in G such that Cj = {vi_i, Vi] for dXXl < i < n. If the two end nodes are the same, i.e., if 
^0 = Vn, the walk is closed. A closed walk is backtrackless if 7^ Cj+i for 1 < i < n—1. A backtrackless 
closed walk is called tailless if e„ 7^ ei. We use the abbreviation TBC for tailless backtrackless closed 
walks. A closed walk with distinct intermediate nodes is called a cycle. The length of the shortest cycle(s) 
in the graph is called girth. In bipartite graphs, including Tanner graphs, all closed walks have even 
lengths, and thus the girth is an even number. 

B. Cyclic Liftings 

Consider the cyclic subgroup of the symmetric group Sn over the set of integer numbers Zjv = 
{0, . . . , A^— 1}, with the group operation defined as composition. The group Gn consists of the circulant 
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permutations defined by — i + d mod N, d e Z^. The permutation tt^ corresponds to d cyclic 

shifts to the right. Corresponding to tt^, we define the permutation matrix J^*^) whose rows are obtained 
by cyclically shifting all the rows of the identity matrix In by d to the left. Clearly, I^^^ = In- There 
is a natural isomorphism between (a) Cn, (b) the group of integers modulo N, Zn, under addition, and 
(c) the group of circulant permutation matrices under multiplication. This isomorphism is defined by the 
correspondence between tt^, d and /^''^ In the following, we refer to d as the shift of the permutation tt^. 

Consider the following construction of a graph G = {V, E) from a graph G = {V, E): We first make 
N copies of G such that for each node v e V^, we have N copies v = . . . , v^~^} in V. For each 
edge e = {u, v} e E, we apply a permutation tt^ e Cn to the N copies of e in £■ such that an edge 
{u^, v^} belongs to £^ iff 7r^(i) — j. The set of these edges is denoted by e. The graph G is called a cyclic 
N -cover or a cyclic N -lifting of G, and G is referred to as the base graph, protograph or projected graph 
corresponding to G. We also call the application of a permutation tt^ e Cjv to the N copies of e, cyclic 
edge swapping, highlighting the fact that the cyclic permutation swaps edges among the N copies of the 
base graph. 

In this work, G is a Tanner graph, and we define the edge permutations from the variable side to 
the check side, i.e., the set of edges e in E corresponding to an edge e = {b, c} G E are defined by 
{b\c^'^^}, i e Zn. Equivalently, e can be described by {h'^'"'''' ^^^\c^}, j e Zn, where (tt*^) ^ is the 
inverse of tt^ in Cn- In other words, if tt^ G Cn is a permutation from variable nodes to check nodes, 
Hd', d' = N — d mod A^, will be the corresponding permutation from check nodes to variable nodes. (Note 
that d' is the additive inverse of d in Zn-) It is important to distinguish between the two cases when we 
compose permutations on a directed walk. 

To the lifted graph G, we associate an LDPC code C, referred to as the lifted code, such that the 
mN X nN parity-check matrix ^ of C is equal to the biadjacency matrix of G- More specifically, H 
consists of a total of mn sub-matrices [II]ij ,1 < i < m, I < j < n, with each sub-matrix of dimension 
N X N, arranged in m rows and n columns. The sub-matrix [II]ij in row i and column j is the circulant 
permutation matrix corresponding to the edge {bj, q} when hij ^ 0; otherwise, [II]ij is the all-zero matrix. 
Let the m X n matrix D — [dij] be defined by [II]ij — 7^*^^ dij e Zn if Kj ^ 0, and dij — -|-oo, otherwise. 
Matrix D, called the matrix of edge permutation shifts, fully describes il and thus the cyclically lifted 
code C- 
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C. ACE Spectrum and ACE Constrained Code Design 

Consider a cycle ^ of length i in a Tanner graph G. The approximate cycle extrinsic message degree 
(ACE) of ^ lfT2]| is defined as ACE(Q = YliA^i '~ 2)' where di is the degree of the ith variable node of ^, 
and the summation is over all the ^/2 variable nodes of ^. Among cycles of the same length, those with 
larger ACE values are less harmful to the performance of iterative decoding algorithms. In [fT2]| . ACE 
constrained LDPC codes were designed whose Tanner graphs were free of short cycles with small ACE 
values. These codes had much better error floors compared to similar random codes. The same idea has 
been used in to devise cyclic lifts for protograph LDPC codes. 

In this paper, we extend the definition of ACE to TBC walks, i.e., for a TBC walk W of length £ in a 
Tanner graph, the ACE is defined as hCE{W) = J^ii'^i ~ 2)^ where di is the degree of the zth variable 
node of W , and the summation is over all the £/2 variable nodes of W . 

Given an LDPC code with a Tanner graph G, the ACE spectrum of depth dmax of G [fT4l is defined 
as a dmax-tuple r]{G) = {r]2, ?72d„,J, where r]2i, i = 1, . . . , dmax, is defined as the minimum ACE 
value of all the cycles of length 2i in G. Note that r]2i = +oo, if there is no cycle of length 2i in G. 
It is desirable to have larger ACE spectrum components for smaller cycles. Given an ACE spectrum 
T] = {ri2, ■ ■ ■ ,V2d„iax) of depth dmax, an LDPC code with a given Tanner graph G (parity-check matrix 
H) is called (dmax, vi) ^CE constrained if for every value of i, 1 <i < dmax, all the cycles of length 2i 
in G have an ACE value larger than or equal to ri2i. In [fT4l| . a generalized ACE constrained construction 
of LDPC codes was introduced, where given an ACE spectrum 77 = (772, •• • iVidmax) of depth dmax, a 
{dmax, v) ACE constrained LDPC code would be designed. Such codes in general outperform the codes 
designed by the original ACE constrained approach of [12] in the error floor region. This is attributed 
to the fact that the codes designed in [fT2l have a "flat" ACE spectrum, while the ACE spectrum for the 
codes constructed in [14] do not have this limitation. The generalized ACE constrained construction was 
also combined with the PEG construction in [fT3l to further improve the error floor. 

Consider an ensemble C"(A, p) of LDPC codes with block length n and degree distribution (A, p). 
A code C from this ensemble is said to have extremal ACE spectrum properties [[T4l if it has an ACE 
spectrum r] = (r]2, . . . , V2dmax), with the property that there is no other code C in C"(A, p) with an ACE 
spectrum rj' = {ri'2, . . . , v'2d„iaJ 7^ ^^'^ — V2i, for 1 < i < dmax- Codes with extremal ACE spectrum 
properties are expected to have good error floor performance [[141 . It is however a formidable task to 
find LDPC codes with extremal ACE spectrum properties or to prove that a code has such properties. 
Nevertheless, it would be desirable to aim at maximizing the ACE spectrum components in the code 
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design process. In this work, we perform this task systematically and in a greedy fashion, focussing on 
one spectrum component at a time, starting from the first ACE spectrum component 772 followed by the 
rest of the components in the increasing order of the cycle lengths until we reach r]2d^ax ■ We refer to the 
proposed scheme as ACE constrained cyclic lifting, or ACE constrained cyclic edge swapping. 



In our construction, we start from a given ACE spectrum 77 = (772, • • • , V2d,nax) of depth dmax, and a 
given base code (graph) C of length n with a certain degree distribution (A, p), and design the non-infinity 
elements of the matrix D of edge permutation shifts such that the lifted code (graph) is {dmax,'n) ACE 
constrained. This results in a {dmax,'n) ACE constrained LDPC code C of length nN, where N is the 
degree of the lifting. The lifted code has the same degree distribution (A, p) as the base code. 

To design a lifting that satisfies an ACE spectrum 77 = (172, • • • , Vid^^ax)^ we need to eliminate all cycles 
of length 2i with ACE values less than r]2i in the lifted graph. To eliminate these cycles in the lifted graph, 
all subgraphs of the base graph whose inverse images include such low ACE cycles should be identified 
and treated. The treatment is performed by the proper selection of the permutation shifts for the edges of 
these problematic subgraphs of the base graph. 

In the following, we show that the problematic subgraphs of the base graph are short TBC walks with 
small ACE values. 

A. Images of Cycles in Cyclic Liftings 

Lemma 1: Consider a cyclic A^-lifting G of a Tanner graph G. Consider a walk W of length ^ in G, 
which starts from a variable node h and ends at a variable node h' with the sequence of edges ei, . . . , e^. 
Corresponding to the edges, we have the sequence of permutation shifts di, . . . ,di. Then the permutation 
shift that maps b, the inverse image of b in G, to b', the inverse image of b' in G, through the walk W is 
d, where 



III. ACE CONSTRAINED CYCLIC LIFTING 




(1) 



i=0 



Proof: The proof is straightforward. See, e.g., [l3]. ■ 

The value of d given in ([T]) is called the permutation shift of the walk from b to b' . Clearly, the 
permutation shift of the walk from b' to b is equal to d' = N — d mod A^. If 6 = b' and all the other 
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nodes are distinct, then the walk will become a cycle and depending on the direction of the cycle, its 
permutation shift will be equal to d or d'. 

Lemma 2: Consider a cyclic A^-lifting G of a Tanner graph G. Consider a TBC walk W of length £ 
and ACE value 7] in G. Suppose that the permutation shift of W is d = 0, and that there is no subgraph 
of W which is also a TBC walk with zero permutation shift. Then, the inverse image ofW'mG consists 
of cycles, each of length i and ACE value r]. 

Proof: Consider a variable node b of the TBC walk W. Consider the copy 6' of 6 in 6 and follow 
the walk in W corresponding to W starting from b\ Since d = 0, the walk ends at b' and is thus closed. 
Since W has no subgraph which is also a TBC walk with zero permutation shift, the above mentioned 
walk never goes through the same node twice. With the same argument, it never meets 6* again prior to 
the end of the walk. The walk is thus a closed walk with distinct intermediate nodes and is therefore a 
cycle. The same argument applies to all the A^ copies of b in b, resulting in A^ non-identical cycles in 
W. The ACE value of each cycle is clearly the same as that of W, as the two have the same variable 
degrees. ■ 

Lemma 3: Any cycle of length i and ACE 7] in the cyclic A^-lifting G of G is projected onto a TBC 
walk W of length i/k and ACE r]/k in G, where /c > 1 is the order of the permutation shifts of W. 

Proof: We first note that the two permutation shifts of W corresponding to the two directions of W 
are the inverse of each other in Z^, and thus have the same order. It is then simple to see that the image 
of a cycle ^ must be a TBC walk W. If the permutation shift of is zero, then k = 1, and the result 
follows from Lemma l2l If 7^ 0, consider a node b of W. Then with an argument similar to the one 
used in the proof of Lemma |2l one can see that the inverse image of W starting from a node b^ E b of ^ 
will return back to 6* after traversing kw edges of W, where w is the length of W. This implies that the 
length i of the resulting cycle is kw. The subgraph ly of G in fact consists of N/k cycles of length £, 
all with the same degree distribution. The ACE result follows readily from the above argument. ■ 
For the cycles, as a special case of TBC walks, we have the following result. 

Lemma 4: Consider a cyclic A^-lifting G of a Tanner graph G. Suppose that ,^ is a cycle of length i 
and ACE 77 in G. The inverse image of ^ in G is then the union of N/ k cycles, each of length ki and 
ACE ki], where k is the order of the permutation shifts of ^. 

Proof: The proof for the number and the length of the cycles in the inverse image of ^ follows from 
the proof for the general case of permutation lifts given in Theorem 2.4.3 of [3J. The ACE results are 
simply a consequence of the fact that each cycle ^' in the inverse image of ^ has exactly the same variable 



degrees as ^ does. The multiplicity of the variable nodes of each degree in ^' however is k times that in 

e. ■ 

In what follows, we refer to the value k in Lemma [3l as the order of the TBC walk, and use the 
notation 0(W) to denote it. It is easy to see that in a cyclic lifting of degree A^, the order of a TBC walk 
W is given by 

^('^' - ^) • P> 

where d is the permutation shift corresponding to W, and gcd denotes the "greatest common divisor." 

Based on Lemmas |2l- SI to eliminate short cycles with low ACE in the lifted graph, we need to make 
sure that short TBC walks with low ACE in the base graph have large orders. 



B. Structure of TBC Walks in the Base Graph 

Consider a base Tanner graph G. Based on the results of the previous section, we are interested in 
enumerating all the short TBC walks in G and making sure that the edge permutation shifts in these 
walks are selected such that the inverse images of the walks are not short cycles with low ACE values. 
The following lemma is simple to prove. 

Lemma 5: A TBC walk in G is either a cycle or consists of at least two cycles. 

Based on Lemma [5l it is clear that the length of a TBC walk in G is at least g. All the TBC walks of 
length g, . . . ,2g — 2 are simple cycles. TBC walks with length > 2g consist of at least two cycles. 

Example 1: Consider the Tanner graph G of Fig. [T] which has 3 variable nodes 61,62,^3 and 3 check 
nodes Ci, 02,03. The graph has two cycles of length 4 and one cycle of length 6 ((7 = 4). All three 
cycles are TBC walks. In addition, there are numerous other TBC walks in G. One example is Wi = 




Cl C2 C3 



Fig. 1. A simple base Tanner graph 

616^026^626^03676365 0264 6263 ciCj^fei, which has length 8. Note that we have used superscripts + and — 
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to denote the direction of the edges in the TBC walk. In the rest of the paper, we may only refer to a walk 
by its sequence of directed edges. For example, Wi can be represented as 65 eje^ef e^l'ejejj'e]^. Other 
examples of TBC walks in G are W2 = e^Cg 65 Cg 6^64 e^e^^, and W3 = 63 Cg 65 64 Cg^Cy Cg^Cg e^^, 
which both have length 10. Note that Wi consists of two cycles of length 4, while each of W2 and W3 
consists of a 4-cycle and a 6-cycle. 

In order to find short TBC walks of the base graph, one can grow a tree from every variable node in 
the graph, as the root, one layer at a time and track the walks on the tree which pass through the root 
node somewhere down in the tree. (Each layer of the tree is constructed by first including all the check 
node neighbors of a variable node except its parent node, and then by including all the variable node 
neighbors of those check nodes except their parent variable nodes.) By the construction, such a walk is a 
backtrackless closed walk. One needs to only select those walks that are also tailless. To find TBC walks 
of length at most 2lmax, one needs to grow the tree up to Imax layers]^ 

C. ACE Constrained Cyclic Edge Swapping Algorithm 

Consider a target ACE spectrum rj{G) = {r]2,r]4, ■ ■ ■ ,V2dmax) for the cyclic A^-lifting G of the base 
graph G. Also, consider a TBC walk W of length w and ACE r](W) in G. The TBC walk W is referred 
to as (potentially) problematic if there exists a divisor k of N such that kw < 2dmax, and kr]{W) < r]kw 
Problematic TBC walks are those whose inverse image in the lifted graph can be cycles that violate the 
target ACE spectrum. One therefore should take proper care in assigning permutation shifts to the edges 
of problematic TBC walks. In particular, if the edge permutation shifts are assigned such that 0{W) = k, 
where k is the positive integer described above, then the kw-th component of the target ACE spectrum 
will be violated by the inverse image of W. The problematic TBC walks can be ordered based on the 
comopnent of the target ACE spectrum that they would violate; the smaller the index of the component, 
the more problematic the walk. In the following, the ordered set of problematic TBC walks in G is denoted 
by W{G). From this set, those that include edge e are denoted by 

Forming an ordered set of problematic TBC walks W(G'), the next step is to go through this set, one 
TBC walk at a time and assign proper permutation shifts to a selected set of edges from the chosen TBC 
walk such that the inverse image of the walk does not violate the target ACE spectrum of the lifted graph. 

^For Tanner graphs with large variable/check degrees, to simplify the search for short TBC walks, one can limit the search to cycles and 
subgraphs that include two cycles. 

'a simpler and still effective approach for ordering the TBC walks in yV(G) is to order them first based on their length, and then for the 
walks of the same length, based on their ACE values. 
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In general, for a problematic TBC walk W, the policy is to select the minimum number of edges of W 
that can make the inverse image of W satisfy the ACE spectrum while maintaining the satisfaction of the 
ACE spectrum by the previously processed problematic TBC walks. 

Example 2: Consider the base graph of Fig. [T] Suppose that the goal is to satisfy the two ACE 
constraints ?74 = ?76 = +00 for the lifted graph. The problematic TBC walks in this case are the 3 
cycles ^1 = e^ejcg e]", ^2 = e^ef Cg ej, and ^3 = e^Cg e^e^Cg . To satisfy the ACE constraint for the 
inverse images of these cycles, we need to make sure that the order of all three cycles is larger than 
one, i.e., the permutation shifts for all three cycles are nonzero. Suppose that the cycles are processed 
in the same order as listed above. To satisfy the ACE constraint for the inverse image of ,^1, we must 
have d{^i) 7^ 0. This can be satisfied by assigning a nonzero permutation shift to only one edge of ^1. 
For example, for = 3, one choice is d{ei) = 1, which results in d(^i) = d{e2) — d{ei) + d{e^) — d{ei) 
mod 3 = 2 7^ 0. Moving on to the inequality 6/(^2) 7^ can also be satisfied by assigning a nonzero 
permutation shift to only one edge of ,^2- To make sure that this assignment will not affect ,^1, we only 
search among the edges of ,^2 that do not belong to ,^1. A proper choice can thus be, e.g., d{e^) = 1, 
which results in d{^2) = 1- Finally, it is easy to see that with the choices for the permutation shifts of ei 
and 65, we have ^(^3) = 1 7^ 0, and thus no more edges need to be processed. 

Given a set of problematic TBC walks W(G) of the base graph G, a target ACE spectrum rj{G) = 
(^2, ^745 • • • 5 V2dmax) for the lifted graph G, and the degree of the cyclic lifting. Algorithm [T] describes 
the proposed ACE constrained cyclic edge swapping. At the output of Algorithm 1, we have the sets 
SwappedSet and ShiftSet, which contain the edges of the base graph that should be swapped, and their 
corresponding permutation shifts, respectively. 

Remark 1: In Algorithm 1, the search for edges to be swapped and the permutation shift assignment 
to these edges are performed in two phases. The first phase is in Steps 4-6, where any edge from 
previously processed TBC walks is removed from the set of candidates for swapping. If the first phase 
fails, in that no edge exists as a candidate for swapping (GandidateSet = 0), then the algorithm switches 
to the second phase in Steps 8-9, where only previously swapped edges are removed from the candidate 
set for swapping. 

Remark 2: To improve the performance of Algorithm 1, in Steps 6 and 9, it is advisable to select edges 
that participate in a larger number of problematic TBC walks. 

Remark 3: By increasing A^, the size of the alphabet space for edge permutation shifts increases. This 
in turn, allows for more problematic TBC walks to be accommodated. In general, one can expect to 
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Algorithm 1 ACE Constrained Cyclic Edge Swapping Algorithm 

Inputs: A target ACE spectrum ri{G) = (+00,774, . . . ■,rj2dmax) of the lifted graph G, an ordered set of 
problematic TBC walks W(G) of the base graph G, and the degree of cyclic lifting A^. 

1) Initialization: ProcessedSet = 0, SwappedSet = 0, ShiftSet = 0. 

2) Select the next problematic TBC walk W G W{G). Denote the length of by w and its ACE value 
by 7]{W). 

3) GurrentSet = edges of W. 

4) GandidateSet = GurrentSet \ ProcessedSet. 

5) If GandidateSet = 0, go to Step 8. 

6) Select the edges £ from the GandidateSet that should be swapped, and assign their permutation shifts 
V from Zn \ {0} such that wO{W) > 2d^ax or 0{W)r){W) > VwOiw)- 

7) SwappedSet = SwappedSet U S, ShiftSet = ShiftSet U V, and ProcessedSet = ProcessedSet U 
GurrentSet. Go to Step 12. 

8) GandidateSet = GurrentSet \ SwappedSet. If GandidateSet = 0, Stop. 

9) Select an edge e from GandidateSet and assign a permutation shift i E \ {0} to it such that the 
inverse images of all the TBC walks in {ProcessedSet) U W satisfy the ACE spectrum ri{G). If this 
is not feasible, go to Step 11. 

10) SwappedSet = SwappedSet U e, ShiftSet = ShiftSet U i. Go to Step 12. 

11) GandidateSet = GandidateSet \ e. If GandidateSet 7^ 0, go to Step 9. Else, stop. 

12) If all the problematic TBC walks in W(G') are processed, stop. Otherwise, go to Step 2. 



achieve a better ACE spectrum as is increased. 

In this work, we use Algorithm 1 in an iterative process to optimize the ACE spectrum of the cyclic 
lifting. The goal is to achieve an ACE spectrum of a certain depth dmax with extremal properties. For 
this, we adopt a greedy approach, where we attempt to maximize the ACE spectrum components, one 
at a time, starting from r]2 and moving towards ri2dmax- Ideally, we would like to achieve ri2i = +00, for 
1 < z < dmax, but this is rarely possible for values of d^ax larger than 3,3 

''One should note that the greedy approach used in this work is not necessarily the best approach in optimizing the ACE spectrum. In 
particular, there may be other achieveable ACE spectrums for a given degree of lifting that result in a better performance for the lifted code. 
Nevertheless, our simulation results demonstrate that the greedy approach of this paper is also quite effective in designing LDPC codes with 
good performance. 
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IV. Numerical Results 

In the following, two examples are presented where the proposed approach is compared with the 
constructions of [[T3l . f[T4l . and [|2l, respectively. 

Example 3: In this example, we consider the construction of a binary irregular LDPC code with rate 
|, and with the degree distribution similar to the one used in the examples of [14]|f . In lfT4l . an LDPC 
code of length n = 1000 was designed using the generalized ACE constrained method. This code had an 
ACE spectrum of (+oo, +oo, 16, 9, 5). 

For a fair comparison, we would like to design an LDPC code with a similar degree distribution and 
with a length of about 1000 by using the proposed method. Since our design is based on the cyclic lifting 
of a protograph, we first need to design a rather small protograph with a degree distribution close to 



that of [[141 rl In this example, we design a protograph with parameters n = 30 and m = 15 by the PEG 
algorithm. The degree distribution of this protograph is X(x) = 0. 2333a;+0. 2250x^ + 0. 1667x^+0. 3750x^^, 
and p{x) = 0.0583a;6 + 0.867x^ + 0.0747a;^ This base graph has ACE spectrum (+oo, 3, 2, 1, 1). 

We then apply the proposed method to design ACE constrained cyclic liftings of various degrees for 
this protograph. The results are shown in Table Hfor lifting degrees = 5, 10, 15, 20, 25, and 30. It can 
be seen from Table |I] that the ACE spectrum of the lifted graph improves by increasing the lifting degree 
N. 

TABLE I 

ACE SPECTRUM OF CYCLIC LIFTINGS OF VARIOUS DEGREES FOR THE (30, 15) BASE CODE OF EXAMPLE[3] 



lifting degree (N) 


ACE Spectrum 


5 


(+00,16,2,2,1) 


10 


(+00,26,2,2,1) 


15 


(+00,26,17,4,2) 


20 


(+oo,+oo,14, 3, 2) 


25 


(+00, +00,17,4,3) 


30 


(+0O, +00,17,9,4) 



The closest block length to 1000 is attained by the lifting degree = 33. The lifted code in this case 

^The degree distribution used in |14|, which is optimized by density evolution, is \{x) = 0.2449a; + 0.20298a::^ +0.00055a;^ +0. 1723a::** + 
0.37923a:*'* and p{x) = x"^ . 

*It is important to note that there is a tradeoff in the selection of the size of the protograph. On one hand, a small size is beneficial in 
having a more compact description of the code and thus a simpler implementation of the encoder and the decoder. On the other hand, it 
may be hard to implement a given optimized degree distribution with high accuracy using a small protograph. A smaller protograph also 
limits the number of variables that are available for the optimization of the ACE spectrum of the lifted code. 
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has length n = 990, and we are able to achieve the ACE spectrum (+00, +00, 17, 10, 5). This improves 
the ACE spectrum of (+cx), +00, 16, 9, 5) obtained in [Ell for n = 1000. 

To have a more fair comparison with the constructions of [TSl and [[T4l|. we also construct two codes with 
length 990 whose degree distributions are the same as our base graph by using the generalized ACE con- 
strained algorithm of [14J and the generalized ACE constrained PEG algorithm of [[T3l FI The constructed 
codes by the methods of [14J and fi3\ have ACE spectrums (+00, +00, 16, 9, 5) and (+00, +00, 17, 9, 4), 
respectively, both inferior to the ACE spectrum of the code designed by the ACE constrained cyclic lifting 
method. Moreover, one should keep in mind that the proposed construction is quasi-cyclic and thus more 
desirable for implementation. 

The frame error rate (FER) and the bit error rate (BER) curves of the three codes over the AWGN 
channel are presented in Fig. |2] The curves are for belief propagation decoding with maximum number 
of iterations 100. As can be seen, the code constructed based on the proposed method outperforms the 
other codes, particularly in the error floor region. 
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— B— GoiicrLili^ed ACE Conslt.iiiieJ PEG [ 13] ; 
























*- 
































V . . . ^ 




1 1 ^1 


1 1.2 1.4 1.6 1.8 2 


2.2 2.4 2.6 2. 







Fig. 2. BER/FER perfoiTnance of the LDPC codes designed in Example [3] 

Example 4: As the second example, we consider the rate-compatible protograph designed in [i2J for 
code rates i, |, |, and |. This protograph is shown in Fig. |3] Different rates are obtained by handling the 
non-trasmitted bits A, B, and C, differently. If all the three bits are constrained to take a zero value, then 
we obtain rate |. This is equivalent to removing the three nodes and their edges from the protograph. 

'More precisely, the codes constructed by methods of 11131 and 1141 have the same variable degree distribution as the base graph, but their 
check degree distributions are slightly different and depend on the details of their design methods. 
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Rate I is obtained if node A is an ordinary non-transmitted node with no bit assignment, but nodes B 
and C are set to zero. Rate | is obtained by setting only node C to zero. For rate |, all three nodes are 
free of any bit assignment. 




Fig. 3. Rate-Compatible protographs of rates 1/2, 5/8, 3/4, and 7/8 used in Example|4l 

To remove the parallel edges, the protograph of Fig. |3]was lifted in [[2l by a lifting of degree 4 using the 
PEG algorithm. This graph, referred to hereafter as the base graph, was then lifted by a cyclic lifting of 
degree 181 using the algorithm of [fTll to obtain an LDPC code of length n = 5792. The ACE spectrums 
of the base graph and the lifted graph of [2J are (+oo, 2, 1, 1, 2) and {+oo, 15, 14, 4, 3), respectively. 

In Table im the ACE spectrum of a number of cyclic liftings of the base graph constructed by the 
proposed method is given. As expected, the spectrum improves as the lifting degree increases. 

TABLE II 

ACE SPECTRUM OF CYCLIC LIFTINGS OF VARIOUS DEGREES DISCUSSED IN EXAMPLE[4] 



lifting degree (N) 


ACE Spectrum 


5 


(+00,15,2,3,2) 


10 


(+CX), 28, 3, 3, 2) 


15 


(+oo,28, 16,2,2) 


20 


(+cxj,+cx),15,3,2) 


25 


(+CX3, +00, 16, 3, 3) 


30 


(+oo,-hoo,16,4,3) 



To compare with the code constructed in [2J, we also design a 181-lifting of the base code using the 
proposed method. The designed graph has the ACE spectrum (-|-oo, +oo, 42, 16, 4), which significantly 
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improves over the ACE spectrum of the code designed in [|2]. We also compare the error rate performance 
of the designed code with that of the code of Q over the AWGN channel in Fig. |4l Belief propagation 
with a maximum number of iterations 100 is used for decoding. Comparing the error rate performances 
of the two codes, one can see that the designed code outperforms the code of [|2|| across the whole range 
of rates, with significantly better performance at the higher rates of 3/4 and 7/8, and particularly in the 
error floor region. 



FER: soUd line, BER: dashed Une 




Fig. 4. BER/FER performance of rate compatible codes designed in Example [4] and those of (2). 



V. CONCLUSION 

In this paper, we propose a method for the construction of finite-length irregular LDPC codes with good 
waterfall and error floor performance. The constructed codes are quasi-cyclic protograph codes and thus 
implementation-friendly. The performance of the codes are enhanced, particularly in the error floor region, 
by the careful selection of edge permutation shifts for the vulnerable subgraphs of the protograph. These 
subgraphs are the ones whose inverse image can be short cycles with low ACE values. We demonstrate 
with a number of examples that the designed codes are superior to previously constructed codes with 
similar parameters in both the ACE spectrum and the error correction performance. 
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